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На основании математической модели технологического процесса показана возможность расчета температуры сырьевой массы 
в барабанных печах. Проведен сравнительный анализ экспериментальных и теоретических данных. 


Неорганические и органические соединения 
фтора имеют большое практическое значение. 
Фтороводород является исходным продуктом для 
получения разнообразных соединений фтора с дру- 
гими веществами. Наибольшее количество безвод- 
ного фтороводорода расходуется в настоящее вре- 
мя на получение элементарного фтора, фтороугле - 
родов, синтетического криолита. 

Промышленный способ производства безвод- 
ного фтороводорода основан на разложении фто- 
рида кальция серной кислотой. Источником СаР 2 в 
процессе разложения является флюоритовый кон- 
центрат, содержащий более 90 мае. % СаР 2 . Рас- 
сматриваемая реакция относится к классу топохи- 
мических эндотермических реакций, поэтому про- 
цесс ведется в обогреваемых вращающихся бара- 
банных печах при 130... 200 °С, куда непрерывно 
подаются флюоритовый концентрат и рабочая 
смесь кислот, состоящая из серной кислоты - 
88. ..96 %, фтороводорода - 1...7 % и воды -2. ..5 %. 

Очевидно, что при управлении технологичес- 
ким процессом необходимо соблюдать жесткие 
требования, предъявляемые как к составу и коли- 
честву сырьевого материала, представляющего со- 
бой гомогенную смесь флюоритового концентрата 
и рабочей смеси кислот, который загружается в ба- 
рабанную печь, так и к тепловым режимам печи. 

Схематическое устройство технологического 
аппарата показано на рис. 1. 

Сырьевой материал поступает в барабанную 
печь через шнековый питатель. Температура в трех 


зонах нагрева барабанной печи поддерживается 
при помощи тэнов, которые могут коммутировать- 
ся как по команде оператора, так и ПЭВМ. 


Реакционный газ 



Для изучения теплового процесса во вращаю- 
щейся барабанной печи и решения вопросов про- 
ектирования автоматизированной системы управ- 
ления технологическим процессом разработана 
статическая тепловая модель, аналогичная описан- 
ной в [1, 6]. В связи с конструктивными особеннос- 
тями технологического аппарата и высокой хими- 
ческой активностью сырья, стационарно располо- 
жить измерительные устройства внутри барабан- 
ной печи не представляется возможным, таким об- 
разом, измерение температуры сырья внутри бара- 
банной печи не производится. Температура сырье- 
вого материала определяется косвенно, по пока- 
заниям датчиков ДТ1, ДТ2, ДТЗ температуры воз- 
духа в соответствующих зонах нагрева. Сложность 
задачи построения математической модели связана 
с низкой степенью корреляции между измеряемы- 
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ми параметрами и показателями качества техноло- 
гического процесса. В статической модели проис- 
ходит вычисление температуры реакционной мас- 
сы в зонах обогрева без учета таких факторов, как 
вращение печи и движение сырья. 

Допущения, принятые при создании математи- 
ческой модели: 

- количество загруженного материала в печь на 
промежутке времени А 1 равно количеству выг- 
руженного; 

- прогрев реакционной массы в зоне обогрева 
происходит равномерно; 

- теплоотдача со вспомогательным технологичес- 
ким оборудованием отсутствует. 
Математический аппарат модели основан на 

уравнениях теплообмена [2, 3, 5]. Условия нагрева 
сырья электрическими тэнами в барабанной печи оп- 
ределяются сложными процессами передачи тепла. 

В силу того, что тэны в рабочем режиме разог- 
реваются до температуры 1000 °С, передача тепла 
излучением в процессе нагрева металлических 
конструкций барабанной печи будет играть реша- 
ющую роль. Расчет теплового излучения основыва- 
ется на законе Стефана-Больцмана. Таким обра- 
зом, лучистая энергия 0„ поверхности тела пло- 
щадью Г определяется из соотношения: 

й,„ = ^ест (1 Г 4 , Вт, 

гдесг 0 -постоянная Стефана-Больцмана, Г-абсолют- 
ная температура черного тела, в - степень черноты. 

Между двумя поверхностями, каждая из кото- 
рых находится в пределах видимости друг от друга, 
происходит непрерывный лучистый теплообмен. 
Результирующий лучистый теплообмен между дву- 
мя излучающими телами зависит от относительной 
интенсивности излучения и коэффициента погло- 
щения каждого из этих тел. Если допустить, что по- 
верхность одного тела серая, а другого тела (или 
окружающей среды) - черная, то тепловой поток 
можно выразить в виде соотношения 

й™ =^ 0 ( 7 І 4 - Г 2 4 ), Вт, (1) 

где Г, и Т г - абсолютная температура излучающего 
и поглощающего тела соответственно. В формулу 
(1) можно ввести поправочные коэффициенты, ко- 
торые учитывают форму тел, углы видимости, а 
также радиационные свойства поверхностей и сре- 
ды, через которую происходит тепловое излучение. 

Учитывая, что скорость вращения барабанной 
печи не превышает 2 об./мин, можно предполо- 
жить, что движение воздуха вокруг печи будет 
иметь ламинарный характер. Таким образом, наг- 
рев металлической конструкции печи осуществля- 
ется также конвективным теплообменом при лами- 
нарном течении. 

При оценке теплообмена между поверхностью 
твердого тела и омывающего ее газа, тепловой по- 
ток <2 К0Ш удобно описывать уравнением: 


а онв = аТАТ, Вт, (2) 

где а - коэффициент теплоотдачи; АТ - перепад 
температуры между поверхностью и газом. 

Выражение (2) описывает закон теплоотдачи 
Ньютона, используемый при анализе всех форм 
конвективного теплообмена. Для ламинарного те- 
чения газа в слое ограниченной толщины вблизи 
поверхности твердого тела теплопередача выража- 
ется через температурный градиент теплоносителя 
в непосредственной близости от поверхности: 
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где А г - коэффициент теплопроводности газа; 
йТ 


сіп 


- температурный градиент в направлении нор- 


мали к поверхности твердого тела. 

Коэффициент теплоотдачи а, называемый также 
коэффициентом теплоотдачи соприкосновением в 
тонком пограничном слое, может быть определён как 
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Коэффициент теплоотдачи сложным образом за- 
висит от геометрии поверхности, характеристик А ж 
теплового потока, а также от физических свойств теп- 
лоносителя. Определение этого коэффициента состав- 
ляет главную проблему конвективного теплообмена. 

Основная задача, которая возникает при реше- 
нии инженерных задач теплопередачи, стоит в том, 
что приходится рассматривать сочетание двух или 
трех основных видов теплообмена. В нашем случае 
тепло поступает от тэнов к барабанной печи как 
конвективным теплообменом, так и излучением. В 
это же время, печь излучает тепловую энергию во 
внешнюю среду и передает тепло реакционной 
массе. С учетом основных положений о теплообме- 
не, изложенных в литературе [2-5] , была построена 
функциональная схема взаимодействия основных 
блоков статической тепловой модели рис. 2. 

Блок БЛ1 предназначен для расчета количества 
тепла, переданного излучением от тэнов к барабану за 
дискретный отсчет времени А/. На входы блока БЛ1 
подаются значения температуры тэнов Г тэн и темпе- 
ратуры Т ст поверхности стали барабана. На основании 
этой информации вычисляется количество тепла, пе- 
реданного излучением от тэнов к барабанной печи. 

В блоке БЛ2 осуществляется расчет количества 
тепла, преданного барабанной печи конвективным 
теплообменом от омывающего её воздуха. На входы 
блока БЛ2 подаются данные о температуре воздуха в 
окрестности камеры, полученные с термопары ДТ1, 
и температуры Т ст поверхности стали барабанной 
печи. На основании этой информации вычисляется 
количество тепла, переданного конвективным теп- 
лообменом. Необходимо отметить, что процесс теп- 
лопередачи может протекать и в обратную сторону, 
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при этом поверхность барабанной печи разогревает- 
ся до температур более высоких, чем температура 
окружающего воздуха, и происходит передача тепла 
от барабанной печи окружающему воздуху. 



Рис. 2. Функциональная схема статической тепловой модели 


В связи с тем, что поверхность барабанной печи 
разогревается до 500 °С, возникает необходимость 
учета тепловой энергии, рассеянной барабанной 
печью во внешнюю среду излучением. Для учета 
этого количества тепла предназначен блок БЛЗ. На 
входы блока БЛЗ подаются значения, полученные с 
термопар о температуре воздуха возле камеры ДТ1 
и температуры поверхности барабана Т ш . На осно- 
вании этой информации вычисляется количество 
тепла, излученного барабанной печь и рассеянного 
во внешнюю среду. 

Блок БЛ4 предназначен для расчета количества 
тепла, переданного от поверхности барабанной пе- 
чи реакционной массе. На входы блока БЛ4 пода- 
ются значения температуры Т ср реакционной массы 
и температуры Т ст поверхности стали барабана. На 
основании этого вычисляется количество тепла, 
переданного реакционной массе излучением и пу- 
тем контактного теплообмена, также в этом блоке 
рассчитывается изменение температуры реакцион- 
ной массы на промежутке времени А 1. 

После определения значений тепловых потоков в 
блоках БЛ1-БЛ4 расчетные значения суммируются. 
Результат суммирования поступает на блок преобра- 
зования, где по суммарному количеству тепла рас- 
считывается изменение ДГ ОТ температуры поверх- 
ности барабанной печи на промежутке времени А/. 

Расчет проводится при условии, что на отрезке 
времени А 1 все параметры технологического про- 
цесса остаются неизменными. 

Для построения математической модели за ос- 
нову была принята функциональная схема, предс- 
тавленная на рис. 2. Период времени А? принят 30 с. 
Модель реализована в табличном процессоре ЕхсеІ. 

На первом этапе исследован процесс прогрева 
барабана в пусковом режиме при отсутствии сырья. 
На рис. 3 показаны результаты моделирования ре- 
жима прогрева технологического аппарата. 



Рис. 3. Режим прогрева технологического аппарата: 1) пока- 
зания датчика ДТ1, 2) расчетные значения температу- 
ры стали барабана 


Расчетные значения температуры поверхности ба- 
рабанной печи к моменту окончания прогрева соотве- 
тствуют регламентной температуре технологического 
аппарата. Таким образом, мы можем перейти к моде- 
лированию рабочего режима барабанной печи, при 
котором сырьё поступает в технологический аппарат. 

На следующем этапе моделирования была пост- 
роена модель с сырьевым материалом. На рис. 4 по- 
казаны сравнительные результаты моделирования. 



Рис. 4. Сравнительные результаты моделирования: 1) тем- 
пература с датчика ДТ1, 2) расчетные значения темпе- 
ратуры поверхности стали барабанной печи, 3) экс- 
периментальные и 4) расчетные значения температу- 
ры сырьевой массы 1-ой зоны 

Рассмотрим более подробно расчетную и экспе- 
риментальную зависимости температуры сырьевой 
массы (рис. 5). 



Рис. 5. Зависимости температуры сырьевой массы в 1-ой зо- 
не барабанной печи: 1) экспериментальные, 2) рас- 
четные значения 
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Рис. 6. Зависимости температуры сырьевой массы в 1-ой зо- 
не барабанной печи: 1) экспериментальные значе- 
ния, пропущенные через фильтр Калмана, 2) расчет- 
ные значения 


СПИСОК ЛИТЕРАТУРЫ 

1. Калинин В.Ф., Матвейкин В. Г., Фролов С. В. Построение пол- 
ной аналитической модели процесса обесфторивания фосфа- 
тов во вращающейся печи. — Тамбов: Изд-во НИИТЭХИМ, 
1987. -45 с. 

2. Уонг X. Основные формулы и данные по теплообмену для инже- 
неров: Пер. с англ. Справочник. — М.: Атомиздат, 1979. —216 с. 

3. Теплообмен излучением: Справочник / Под ред. А.Г. Блоха, 
Ю.А. Журавлева, Л.Н. Рыжкова. — М.: Энергоатомиздат, 1991. 
-432 с. 


В целях исключения нежелательных выбросов, 
связанных с методикой проведения эксперимента, 
экспериментальные данные были пропущены через 
фильтр Калмана. Сглаженная экспериментальная и 
расчетная зависимости представлены на рис. 6. 

Экспериментальные и расчетные зависимости 
имеют расхождение не более 1 %. 

Сравнительный анализ расчетных и экспери- 
ментальных данных (рис. 6) позволяет заключить о 
достаточно адекватной математической модели 
теплового процесса барабанной печи. Всё это дела- 
ет возможным использовать разработанную модель 
теплового процесса барабанной печи для управле- 
ния технологическим процессом получения без- 
водного фтороводорода. 

4. Кудинов А.А., Кудинов В.А. Теплообмен в многослойных 
конструкциях. Инженерные методы. — Саратов: Изд-во Сара- 
тов. ун-та, 1992. — 136 с. 

5. Справочник по теплообменным аппаратам / Под ред. П.И. Ба- 
жана, Г.Е. Каневец, В.М. Селивестрова. — М.: Машинострое- 
ние, 1989. —329 с. 

6. Арутюнов В.А., Бухмиров В. В., Крупенников С.А. Математи- 
ческое моделирование тепловой работы промышленных пе- 
чей. — М.: Металлургия, 1990. —229 с. 


ѵпк тел 1 -Л Ітч яі 





